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Abstract 

Recent  theoretical  work  by  Hall  &  Seddougui  (1989)  has  shown  that  strongly  non¬ 
linear,  high  wavenumber  Gortler  vortices  developing  within  a  boundary  layer  flow  are 
susceptible  to  a  secondary  instability  which  takes  the  form  of  travelling  waves  confined 
to  a  thin  region  centred  at  the  outer  edge  of  the  vortex.  This  work  considered  the 
case  in  which  the  secondary  mode  could  be  satisfactorily  described  by  a  linear  stability 
theory  and  in  the  current  paper  our  objective  is  to  extend  this  investigation  of  Hall 
&  Seddougui  (1989)  into  the  nonlinear  regime.  We  find  that  at  this  stage  not  only 
does  the  secondary  mode  become  nonlinear  but  it  also  interacts  with  itself  so  as  to 
modify  the  governing  equations  for  the  primary  Gortler  vortex.  In  this  case  then,  the 
vortex  and  the  travelling  wave  drive  each  other  and,  indeed,  the  whole  flow  structure  is 
described  by  an  infinite  set  of  coupled,  nonlinear  differential  equations.  We  undertake 
a  Stuart- Watson  type  of  weakly  nonlinear  analysis  of  these  equations  and  conclude,  in 
particular,  that  on  this  basis  there  exist  stable  flow  configurations  in  which  the  trav¬ 
elling  mode  is  of  finite  amplitude.  Implications  of  our  findings  for  practical  situations 
are  discussed  and  it  is  shown  that  the  theoretical  conclusions  drawn  here  are  in  good 
qualitative  agreement  with  available  experimental  observations. 
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§1.  Introduction. 


Our  concern  is  with  an  asymptotic  description  of  the  three-dimensional  break¬ 
down  of  steady,  spanwise  periodic,  large  Gortler  vortices  via  nonlinear  travelling  waves 
and  is  a  natural  extension  of  the  work  of  Hall  &  Seddougui  (1989)  who  investigated 
the  current  problem  in  the  context  of  breakdown  through  very  small  waves. 

Experiments  by  Bippes  &  Gortler  (1972)  and  Aihara  &  Koyama  (1981)  suggest 
that  the  breakdown  of  vortices  leads  to  a  flow  with  wavy  boundaries  and  which  is 
periodic  in  time.  Later,  Koliama  (1987)  studied  vortex  instabilities  in  boundary  layers 
over  a  laminar  flow  control  wing,  and  observed  that  a  secondary  instability  of  Gortler 
vortices  appeared  which  was  confined  to  a  thin  zone  at  the  top  of  the  extent  of  vor¬ 
tex  activity.  Additionally,  he  noted  that  the  instability  propagated  with  a  speed  that 
approached  the  free  stream  speed  as  the  disturbance  travelled  downstream,  and  the 
onset  of  the  unsteady  flow  perturbations  was  eventually  followed  by  transition  to  tur¬ 
bulence.  Peerhossaini  &  Wesfreid  (1988)  studied  the  fluid  flow  in  a  boundary  layer 
in  a  concave  section  of  a  curved  channel.  Like  Kohama,  they  found  that  when  the 
secondary  breakdown  of  Gortler  vortices  appeared  it  did  so  at  the  top  of  the  vortices. 
Further  downstream  they  reported  another  instability,  but  which  was  confined  to  a 
zone  close  to  the  wall. 

Much  of  the  early  theoretical  work  concerned  with  Gortler  vortices  concentrated 
on  the  linear  stability  of  external  flows  over  concave  walls.  Notable  contributions  were 
made  by,  among  others,  Gortler  (1940),  Smith  (1955),  Hammerlin  (1956)  and  Floryan 
&  Saric  (1979).  However,  Hall  (1982a,  b,  1983)  demonstrated  that  much  of  this  early 
work  was  fundamentally  flawed  for  all  the  previous  analyses  had  invoked  the  parallel 
flow  approximation  (which  essentially  assumes  that  the  basic  flow  in  which  the  vortices 
lie  is  independent  of  the  streamwise  co-ordinate  and  so  neglects  the  effect  of  bound¬ 
ary  layer  growth).  This  approximation  enables  the  linear  instability  equations  to  be 
expressed  as  ordinary  differential  equations  but  Hall  illustrated  that  this  assumption 
is  unjustifiable  except  in  the  limit  of  smaL  vortex  wavelength,  and  indeed  this  is  the 
explanation  for  the  considerable  inconsistencies  in  the  results  of  the  previous  studies. 
Additionally,  in  the  case  of  small  vortex  wavelength  the  Gortler  instability  may  be 
described  by  an  asymptotic  structure  which  takes  account  of  boundary  layer  growth 
in  a  fully  rational  manner  and  thence  the  parallel  flow  approximation  is  rendered  su¬ 
perfluous  in  the  only  situation  in  which  it  has  any  relevance  whatsoever.  For  larger 
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wavelengths  the  problem  is  fully  non-parallel  and  the  linear  stability  equations,  which 
now  take  the  form  of  a  set  of  partial  differential  equations,  have  to  be  solved  numeri¬ 
cally,  see  Hall  (1983).  This  paper  showed  two  significant  features  of  this  non-parallel 
flow  problem,  namely  that  the  ideas  of  a  unique  neutral  stability  curve  and  of  unique 
growth  rates  at  a  specified  downstream  location  are  inapplicable  to  the  Gortler  prob¬ 
lem  because  the  location  where  a  vortex  commences  to  grow  is  dependent  upon  the 
position  and  the  size  and  shape  of  the  imposed  disturbance. 

Questions  concerning  the  development  of  nonlinear  non-parallel  vortices  within 
growing  boundary  layers  were  answered  by  Hall  (1988).  This  numerical  investigation 
showed  that  as  the  nonlinear  disturbance  evolved,  the  perturbation  energy  became 
concentrated  in  the  fundamental  and  mean  flow  correction;  a  conclusion  consistent 
with  the  weakly  nonlinear  theory  of  Hall  (19826)  valid  for  small  wavelength  vortices.  It 
is  well  known  that  Gortler  vortices  set  up  in  an  experiment  conserve  their  wavelength 
as  they  move  downstream.  Since  the  boundary  layer  itself  thickens  it  follows  that 
the  local  nondimensional  vortex  wavenumber  becomes  large  as  the  vortex  develops. 
Thus  the  small  wavelength  limit  in  the  external  Gortler  problem  is  appropriate  to 
the  ultimate  development  of  any  fixed  wavelength  vortex,  and  hence  sufficiently  far 
downstream  in  many  flows  the  asymptotic  work  of  Hall  (1982a,  6)  becomes  applicable. 

As  with  all  weakly  nonlinear  investigations,  the  results  of  Hall  (19826)  are  valid 
only  within  the  neighbourhood  of  the  point  where  the  considered  perturbation  is  neu¬ 
trally  stable.  For  vortices  of  wavenumber  e-1,  0  <  e  «  1,  their  development  down¬ 
stream  of  the  point  of  neutral  stability  is  dictated  by  the  solution  of  a  pair  of  coupled 
nonlinear  partial  differential  equations  which  adopt  a  simple  asymptotic  structure  at 
large  values  of  X,  where  eX  {X=  0(1))  denotes  the  distance  of  the  vortex  down¬ 
stream  of  the  neutral  point.  Formally,  for  large  X,  Hall  Sc  Lakin  (1988)  showed  that 
this  asymptotic  structure  could  be  used  to  deduce  the  flow  configuration  for  fully 
nonlinear  Gortler  vortices,  at  which  point  the  mean  flow  correction  generated  by  the 
presence  of  the  vortices  becomes  as  large  as  the  basic  (undisturbed)  flow  itself.  These 
fully  nonlinear  vortices,  as  described  by  Hall  Sc  Lakin  (hereafter  referred  to  as  HL) 
are  of  the  type  whose  stability  to  travelling  wave  disturbances  will  be  considered  here. 
For  this  reason,  it  is  convenient  to  now  briefly  describe  the  main  results  of  HL. 

Consider  a  vortex  of  wavelength  27re  developing  in  a  boundary  layer  flow  at  a 
Gortler  number  0(€-4),  and  suppose  that  the  flow  is  neutrally  stable  ai  the  down- 


stream  position  xn.  For  x  >  xn,  HL  demonstrated  that  the  flow  field  assumes  the 
structure  illustrated  in  Figure  1,  where  the  vortex  activity  is  restricted  to  region  I  and 
decays  exponentially  within  thin  shear  layers  IIa,b.  Suppose  that  these  are  located  at 
7/i  and  7/2 -  Then  exterior  to  the  shear  layers,  in  zones  IIIa,b ,  the  vortex  is  annihilated 
and  the  mean  flow  satisfies  the  usual  boundary  layer  equations.  In  the  central  part  of 
the  flow,  region  I,  the  mean  flow  is  dictated  by  a  solvability  condition  on  the  funda¬ 
mental  component  of  the  vortex  and  adjusts  itself  so  as  to  make  the  fundamental  and 
all  harmonics  finite  within  I.  The  vortex  is  determined  by  the  mean  flow  equations, 
and  the  positions  of  the  shear  layers  y\,V2  change  as  they  move  downstream:  in  fact 
they  are  given  by  a  free  boundary  value  problem  derived  from  the  boundary  layer 
equations.  Recall  that  this  small  wavelength  asymptotic  approach  is  relevant  to  many 
flows  of  practical  importance  for  the  reasons  given  previously. 

We  shall  study  secondary  instabilities  of  the  flow  in  the  shear  layers  IIa,b.  Per¬ 
turbations  which  are  spanwise  periodic  travelling  waves  are  imposed  on  the  flow  within 
these  zones  and  their  development  monitored.  These  disturbances  are  tt/2  radians  out 
of  phase  with  the  fundamental  in  the  spanwise  direction,  so  that  any  secondary  insta¬ 
bilities  which  occur  produce  locally  wavy  vortex  boundaries  in  IIa,b.  At  the  outset 
it  is  not  clear  that  the  shear  zones  should  be  susceptible  to  the  wavy  modes  but  Hall 
&  Seddougui  (1989)  (hereafter  referred  to  as  HS)  presented  an  argument  as  to  why 
this  should  be  the  cas„e.  Basically,  they  considered  a  model  equation  which  governs 
the  nonlinear  growth  of  time-dependent  Gortler  vortices  in  curved  channel  flows.  The 
instability  of  a  steady,  streamwise  independent  vortex  to  travelling  vortex-like  pertur¬ 
bations  was  investigated  and  it  was  demonstrated  that  if  the  travelling  perturbations 
were  in  phase  with  the  fundamental  component  of  the  steady  vortex  then  the  flow  is 
always  stable  whereas  if  the  disturbance  were  7r/2  out  of  phase  then  instability  could 
result.  Moreover,  if  the  Gortler  number  was  large,  the  disturbance  was  confined  to 
the  shear  layers  and  decayed  away  from  the  centres  of  these  layers.  This  analogy  with 
the  curved  channel  flow  problem  led  HS  to  analyse  disturbances  in  the  boundary  layer 
problem  which  were  concentrated  in  the  shear  layers  I  la,  b.  As  the  disturbances  de¬ 
cayed  exponentially  away  from  the  centre  of  these  shear  layers,  HS  deduced  that  unless 
exponentially  small  terms  were  to  be  matched  the  two  shear  layers  could  be  treated 
independently. 

In  HS  the  wavy  perturbations  were  assumed  to  have  an  infinitesimally  small  am- 
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plitude.  It  was  found  that  in  this  case  these  wavy  disturbances  are  fixed  by  a  pair 
of  coupled,  second  order,  linear  ordinary  differential  equations  which  were  solved  nu¬ 
merically.  Hall  &  Seddougui  found  only  one  neutrally  stable  travelling  wave  solution 
of  these  equations  and  this  information  was  used  to  explain  some  of  the  experimental 
observations  of  wavy  vortices;  notably  those  of  Koliama  (1987)  and  Peerliossaini  & 
Wesfreid  (1988).  However,  later  studies  by  Bassom  &  Seddougui  (1989)  modified  the 
conclusions  of  HS  by  demonstrating  that  infinitely  many  neutrally  stable  modes  are 
possible  and  properties  of  these  modes  in  a  high-wavenumber  limit  were  described 
asymptotically. 

Here  we  aim  to  increase  the  amplitude  of  the  wavy  disturbances  to  extend  the 
ideas  of  HS  to  the  nonlinear  regime.  Importantly,  it  is  found  that  at  the  point  at 
which  the  wavy  mode  becomes  nonlinear,  the  steady  vortex  flow  as  determined  by 
HL  is  affected  by  self-interactions  of  the  wavy  mode.  Additionally,  at  this  stage 
all  the  harmonics  of  the  wavy  mode  become  sufficiently  large  so  as  to  influence  the 
fundamental  at  leading  order.  Consequently,  we  found  that  this  nonlinear  problem  is 
governed  by  a  fully  coupled,  infinite  set  of  differential  equations  which  necessarily  have 
to  be  solved  numerically.  Rather  than  tackle  this  complete  system,  we  concentrate  on 
a  weakly  nonlinear  limit  which  enables  further  analytical  progress  to  be  made.  The 
results  of  this  work  are  considered  in  the  light  of  the  available  experimental  evidence 
and  conjectures  made  concerning  the  results  which  could  be  anticipated  from  a  strongly 
nonlinear  calculation. 

The  procedure  adopted  in  the  remainder  of  the  paper  is  as  follows.  In  section  2 
we  formulate  the  problem  and  develop  the  defining  equations  for  the  nonlinear  wavy 
modes  in  section  3.  We  analyse  the  small  amplitude  limit  of  these  equations  in  section 
4  and  in  section  5  we  discuss  our  results,  consider  the  full  numerical  solution  of  the 
problem  and  draw  some  brief  conclusions. 


§2.  Formulation  of  the  problem. 

The  flow  considered  here  is  that  as  described  in  HL.  We  study  the  flow  of  an 
incompressible,  viscous  fluid  of  density  p  and  kinematic  viscosity  v  over  a  wall  of 
varying  curvature  ~x  {X/L).  Here  R  is  a  typical  radius  of  curvature,  X  denotes 
distance  along  the  wall  and  L  is  a  typical  lengthscale  along  the  wall.  Much  of  the 
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following  detail  concerning  tlie  basic  flow  structure  in  the  absence  of  the  imposed 
travelling  disturbances  follows  very  closely  that  described  in  HL  and  HS  and  so  we 
keep  our  description  of  this  aspect  as  brief  as  possible  consistent  with  the  aim  of 
retaining  clarity. 

The  Reynolds  number  Re  for  the  flow  is  defined  by 


(2.1) 


where  Uq  is  a  typical  fluid  velocity.  The  curvature  parameter  6  is  defined  by 


6  = 


L 

72’ 


(2.2) 


and  our  concern  is  with  the  case  in  which  Re  — >  oo  whilst  the  Gortler  number  G, 
given  by 

G  =  2  Re*  6,  (2.3) 

maintains  a  constant  0(1)  value.  We  denote  time  by  T  and  co-ordinates  along  the  ■wall, 
normal  to  the  wall  and  in  the  spanwise  direction  by  X,  Y  and  Z  respectively.  If  the 
corresponding  velocity  vector  is  (U,  V,  W)  then  we  choose  nondimensional  co-ordinates 
(x,y,z)  and  velocities  (u,v,w)  given  by 


{x,y,z)  =  L~l  (x,Re'’Y,Re?z)  , 


( u,v,w )  =  Uq 1  (u,Re’V,Re>w)  .  (2.4a, 6) 


In  the  work  considered  here  we  restrict  ourselves  to  examining  flows  for  which 
u  — ►  1  as  y  — >  oo,  although  there  is  no  formal  difficulty  in  extending  our  treatment 
to  account  for  more  general  flows.  We  write  the  pressure 

P=§P,  (2.4C) 


and  substitute  the  expressions  (2.4)  into  the  governing  continuity  and  Navier-Stokes 
equations.  Defining  the  nondimensional  time  variable  by  t,  where  t  =  UqL'^T,  we 


find  that  when  terms  of  order  Re~  *  are  neglected  in  the  governing  equations,  these 
equations  assume  the  forms 

du  dv  dw  . 

fc  +  dy  +  T)z=0'  (2-5a) 

du  du  du  du  d2u  d2u  . 

Tt+ateJrVaZ+wTz  =  W+eZ'  (  ) 
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dv 


dv 


dv 


dv 


^7  +  uir:  +  v—  +  w—  =  - -GXu *  - 


dp  d2v  d2v 


dt 


dx 


dy 


dz 


+ 


+ 


dy  dy*  dz2' 


dw  dw  dw  dw  dp  d2w  d2w 

sT  +  ”  to  +  voj  +  ”  al  =  ~Tz  +  W  +  ^ 


(2.5  c) 


(2.5d) 


The  motivation  for  the  work  contained  in  HL  was  to  obtain  a  solution  of  (2.5) 
which  satisfies 


u  =  v  =  w  =  0  on  y  =  0,  u  — >  1  as  y  — »  oo,  (2.6) 

and  which  is  an  asymptotic  solution  valid  in  the  limit  of  small  vortex  wavelength.  The 
Gortler  number  G  expands  according  to 

G  =  Gq(~  4  +  G\e~  3  +  . . . ,  (2.7) 

where  e~ 1  (0  <  e  <C  1)  is  the  wavenumber  of  the  vortices.  In  the  core  part  of  the 
boundary  layer  (region  I),  the  leading  order  ^-independent  part  of  the  downstream 
velocity  component  u(x,y)  satisfies 

du 

°0XU—  =  1,  (2.8) 

and  this  flow  is  driven  by  the  finite-sized  steady  vortices.  These  vortices  produce  a 
mean  flow  correction  term  which  is  as  large  as  the  basic  flow  which  would  exist  in  the 
absence  of  the  vortices  and  so  the  total  mean  flow  field  is  altered  at  leading  order  by 
the  presence  of  the  vortices.  The  calculation  for  the  vortices  within  the  core  region 
(see  Figure  1)  reveals  that  the  vortices  are  trapped  between  y\  and  t/2  (the  values  of 
which  are  found  in  the  course  of  the  computation)  and  decay  to  zero  within  the  thin 
shear  layers  I  la,  b. 

The  work  given  by  Davey  et.  al.  (1968)  showed  that  Taylor  vortex  flow  is  unstable 
to  perturbations  differing  in  phase  from  the  fundamental  component  of  the  steady 
vortex  flow  by  90°.  In  HS,  the  authors  sought  a  secondary  instability  that,  produced 
locally  wavy  vortex  boundaries  in  the  shear  layers  JIa,  b ;  akin  to  the  wavy  surfaces 
predicted  by  Davey  et.  al.  in  the  case  of  the  secondary  instability  of  their  Taylor  vortex 
flow.  In  keeping  with  this  aim,  HS  examined  solutions  for  out  of  phase  perturbations 
proportional  to 

E  =  cxj>(^^J  K (x)dx  —  >  (2.9a) 
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where  the  wavenumber  K  was  expanded  as 

K  =  K0  +  e$K1  +  ...,  (2.96) 

and  fl  is  the  constant  frequency.  The  wavy  vortex  mode  was  assumed  to  be  a  short 
wavelength,  high  frequency  mode  which  was  necessary  in  order  to  ensure  that  this 
secondary  instability  was  to  remain  trapped  within  Ila  or  lib,  an  assumption  which 
is  consistent  with  the  available  experimental  observations.  The  scales  in  (2.9)  are 
chosen,  from  Hall  (19826),  so  that  ~  and  |j-  +  -  0  in  the  shear  layers.  The 

latter  scaling  ensures  that  the  waves  travel  with  the  speed  of  the  fluid  in  the  shear 
layer,  as  observed  in  the  work  of  Kohama  (1987). 

In  HS  the  secondary  instabilities  were  supposed  to  be  of  infinitesimal  amplitude 
so  that  the  steady  vortex  flow  was  unaffected  by  the  presence  of  these  wavy  modes. 
Our  present  aim  is  to  allow  this  travelling  vortex  to  be  of  sufficiently  large  amplitude 
so  that  self-interactions  of  this  mode  do  modify  the  governing  equations  for  the  steady 
vortex.  We  then  have  the  case  of  a  fully  coupled  problem  in  which  the  primary  vortex 
and  the  wavy  mode  must  be  determined  simultaneously;  in  contrast  to  the  calculation 
presented  in  HS  in  which  the  wavy  mode  was  found  only  after  the  steady  vortex  had 
been  fully  determined.  In  the  following  section  we  develop  the  asymptotic  expansions 
for  both  the  vortex  and  the  nonlinear  wavy  mode  under  the  assumption  that  the  latter 
is  confined  to  within  the  shear  layer  Ila,  an  analysis  which  may  be  trivially  modified 
to  account  for  modes  wit!  'a  the  other  shear  layer  if  so  desired. 


§3.  The  asymptotic  structure  of  the  wavy  modes. 


In  HL  it  was  demonstrated  that  the  layers  IIa,b  are  of  thicknesses  0(e*)  and  so 
in  Ila  we  define  the  co-ordinate 

( V  ~  V2 ) 


£  = 


e* 


(3.1) 


where  y2{x)  is  the  location  of  the  zone  Ila.  Hence,  in  this  layer  we  replace  d/dx  and 
djdy  by 

_  J/2  ,  J__#_ 

dx  dt  el-  dt 

respectively.  The  vortex  activity  is  wholly  within  the  core  region  I  between  Ila  and 
lib,  and  the  algebraically  decaying  vortices  in  I  are  reduced  to  zero  exponentially 
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within  the  shear  layers.  In  the  analysis  described  below  we  find  that  the  leading  order 
mean  flow  terms  across  the  whole  of  the  boundary  layer  are  unchanged  by  the  presence 
of  the  nonlinear  wavy  mode  terms.  Additionally,  the  steady  vortex  terms  within  the 
core  region  I  are  unaffected  by  the  wavy  mode,  but  this  latter  disturbance  does  alter 
the  vortex  quantities  within  the  shear  layers. 

We  consider  the  amplitude  of  the  wavy  mode  to  be  such  that  the  equations  given 
in  HS  to  determine  the  steady  vortex  within  the  shear  layers  are  modified  due  to 
interactions  of  the  wavy  disturbances.  It  is  found  that  when  the  wavy  mode  is  as  large 
as  the  steady  vortex  (=  0(e>),  by  HL)  this  change  in  the  character  of  the  governing 
equations  occurs  and  we  are  led  to  consider  total  velocity  and  pressure  fields  within 
the  shear  layer  Ila  which  take  the  forms 


u  =5q  +  §1  +  eU2  +  e*u3  +  ...  +  [c^beO  +  e* CUoei  +  . . .] 


+ 


4 

+  e* 


E  (Cuocn  +  Suo,n)En  +  e*  E  (Cuicn  +  Su\,n)En  +  . . .  +  c.c. 

n=l  n=l 

oo 

E  U On  En  +  . . .  +  c.c.  +  .  .  .  , 

n— 1 


(3.2a) 


v  =0q  +  Si  +  . . .  +  c  »  [cVqcO  +  e*  CVq ci  +  . . .] 


+  e“» 


E  ( Cv0cn  +  Sv0tn)  En  +  E  {CV1  en  +  Sv  1<n)  E*  +  . , 

n=l  n=l 


.  +  C.C. 


+ 


E  V On  En  +  . . .  +  c.c. 

n=l 


+  . . .  , 


(3.26) 


V) 


=e~  1  [sWo.-O  +  er5Wb.*l  +  •  • .] 


,  _  t 

+  e  * 


E  (C'WOen  +  Swo,n)En  +  e*  E  (Cwien  +  Sw\in)  En  +  . . .  +  c.c. 

“n=l  n=l 


E  W'O nEn  +  ...  +  C.C. 

n—  1 


■f  • . . , 


(3.2  c) 
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=e  4Po  +  e~  *3°  Pi  +  •  •  •  +  e”  *  [ePocO  +  c»  CPoci  +  . . .] 


oo  oo 

+  Yl  (Cpocn  +  Spo,n)En  +  X  {Cpien  +  Spi,n)En  +  ...  +  c.c. 

-  _ -»  _ 1 


+  £  3  P0nEn  +  ♦  . .  +  C.C. j  +  .  .  ♦  . 

(3.2d) 

In  (3.2)  C  and  5  denote  cos  (2/e)  and  sin(^/e)  respectively,  all  unknown  coeffi¬ 
cients  are  functions  of  x  and  f  and  c.c.  denotes  complex  conjugate.  The  expansions 
are  more  complicated  versions  of  those  presented  in  HS,  and  we  recall  that  the  vari¬ 
ation  E  defined  by  (2.S)  dictates  the  streamwise  development  of  those  parts  of  the 
total  velocity  and  pressure  fields  corresponding  to  the  wavy  mode.  As  is  expected, 
when  the  amplitude  of  the  wavy  mode  is  decreased  from  that  value  implied  in  (3.2), 
we  retrieve  the  expansions  of  HS.  We  remark  that  higher  harmonics  of  the  spanwise 
variation  lead  to  terms  of  smaller  sizes  than  those  given  in  (3.2),  see  HL.  These  terms 
do  not  enter  the  determination  of  the  primary  vortex  or  the  wavy  mode  and  so  have 
been  omitted  in  (3.2).  However,  we  do  observe  that  all  harmonics  of  the  spatial  and 
temporal  variation  E  must  be  accounted  for  since  they  all  play  nontrivial  roles  in  the 
subsequent  analysis. 

We  find  the  unknowns  in  (3.2)  by  substituting  these  forms  into  (2.5)  and  equating 
coefficients  of  specific  terms.  The  leading  order  mean  flow  quantities  So  and  Si  are 
found  to  be  precisely  as  in  HL  so  that 


.  /  a  +  2ijz 

~  V  1 


V  <?ox(a  +  21/2) 


Here  a(x)  is  an  arbitrary  function  which  may  be  determined  by  considering  the  so¬ 
lutions  in  region  I,  see  HL.  In  HL,  the  fundamental  vortex  term  Vqcq  was  found  by 
deriving  a  solvability  criterion  on  higher  order  vortex  t  5  within  the  shear  layer 
Ila.  With  the  large  amplitude  wavy  disturbance  imposed,  this  solvability  condition 
assumes  the  form 


3^  +  23-  +  f = 

a  +  27/2 

oo  00  f 

X  (w'on  (V0sn  -  GoxSoWO.n)*  +  C,CV  “  ^OX^O X  (  t>0o» 
«  _ 1  ' 


5  *  \ 

^♦4 


(3.4a) 
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where 


du3  ,,  , 

aT  = /(I) 


x'(a  +  2t/2) 


+ 


£- 


Vn2. 


OcO 


+ 


!/(£( 

n=l 


3G0X2  >/Cox(«  +  22/2)''  2y/ Gox(a  +  27/2) 

dzi*  dum  \  \ 

VOcn  +  UO.n  ■  +  WQcnUOin  ~  ^O.nUQen  +  C.C.J  J  dZ‘ 


(3.46) 

These  equations  are  derived  following  procedures  identical  to  those  described  in 
HS.  Here  b(x)  is  a  function  arising  from  the  solution  in  the  core,  f(x)  is  a  function 
which  can  only  be  found  by  considering  higher  order  terms  (which  proves  unnecessary 
here)  and  an  asterisk  on  a  quantity  denotes  the  complex  conjugate  of  that  quantity. 
The  equations  equivalent  to  (3.4)  for  the  case  of  infinitesimally  sized  wavy  vortices, 
as  in  HS,  may  be  retrieved  by  removing  the  terms  comprising  products  of  wavy  mode 
perturbation  quantities  in  (3.4a,  6). 

Comparing  coefficients  of  sin  (2/e)  in  (2.5 a  —  c)  and  of  cos  (2/e)  in  (2.5d)  yields 

the  following  relations  between  the  leading  order  wavy  vortex  terms.  We  obtain 

Ov  Qln 


in  (KqUq  -  Sl)uo,n  +  =  ~«(W, 


dUi 

w 


—  wo cn  =  0,  (3.5a) 

(3.56) 
(3.5c) 
(3.5  d) 

(3.5e) 

,  «0en  = - Qf~V°cn,  POcn  = - (3.5/) 

Substituting  into  (3.4)  and  simplifying  reveals  that  the  steady  fundamental  vortex 
term  VbcO  satisfies 


in  ( A'o«o  -  SI)  vo,n  +  GoxHouo.n  =  -vo.n, 

in  ( KoUo  -  SI)  WQen  =  -POsn  ~  W0cn. 
Now  (3.56,  c)  are  only  consistent  if  Ktf&O  =  ^  and  then  we  have 


WOcn  = 


dvo.n  dui 

^0*n  —  ^Omj 


d$  ’  “v‘"  d£ 

By  an  entirely  similar  procedure  we  can  obtain 

3U0cn  55l 

WQ,n  = 


POtn  — 


dvo,n 
■  di  • 

“W0cn 


32Vh.n  f  y’  (a  +  b 

-r(yMiSl  — 


QxVG0X  3  3(a  +  2t/2)J 


vl 


0  iO 


^(e(N.»I2  +  K.I2))- 

n=l 


1 

6  n=l 


fV  G0x(a  +  2y2) 


Qrf 

2WonuO*n  -  y/ Gox  (a  +  2y2)~d£-V0cn  +  C.C. 


VbcO 


(3.6) 
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Examining  terms  proportional  to  En  in  (2.5 d)  vields  the  governing  equation  for 
the  mean  flow  correction  term  W o„ ; 


d2W 


On 


_  Clui 
—  tn  |  K\uq  +  — — 
uo  J 

d2vo,n  d2V0, 


W  0„  = 


1  ( r  r  ^"^O.n  0~VQeQ 

2  ^*0e0  nr 9  V0sn 


di2 


dt2 


\  ,1  V'  d  (  ^O.(n-k)  5t>0e(n-k)^ 

/  2  ^  \  0e*  <9£  VOak  d£  /  ' 

Jb=-  oo 


(3.T) 

To  obtain  the  equation  for  the  wavy  vortex  term  VQsn  we  need  to  consider  the 
governing  equations  for  the  higher  order  terms  in  the  scheme  outlined  in  obtaining 
(3.5).  Following  this  procedure  we  find  from  (2.56,  c)  that  ui4n  and  t>i4n  satisfy 


dus  2  C^UO,n 

in  (A'ofll  +  /fl«o)«0*n  +  VOin*^-  -  W()»^OcO - 

,  V1  f  dZ7()(n-Jt)  TJr  ^ 

+  2s  ^ - -  W^OfcttOe^-fcjj  =  “«1*»  -Vim—, 

Jk=-  oo 


(3.8a) 


and 


in  {KqUI  +  Kiuo)  VQtn  -  “  ^beOWOn  +  G()X“1U0»» 

oo 

-  2j  vockWo(n_k>)  =  -GoxUoui,n  -  Vi,n. 

i fe=—  OO 


(3.86) 


Now  using  the  results  (3.5 e,/)  we  find  that  (3.8a, 6)  are  consistent  for  the  unknowns 
(«l.n,«i.n)  only  if 


d2V\ 


o.» 


+  £ 


1  _  x'  (fl  +  2 y2) r  _  6 

0£2  ,sL3(a  +  2y2)  9Xv/^X  3 


*>0»n 


2m 


ATl-s/a  +  2j/2 


^0»n 

U0,n  =  - - 


Vn2. 


OcO 


3  *•  y/m  a  4-  2y2 

+  yp  £  (i^'Ocni2  +  j^Oinl2)  -  C?OX(c  +  2T/2)uo,n  “ 
n=l 


/. 

3 

Q^O(n-k) 


(3.9a) 


-  j  ]C  1  v/Gox(a  +  2j/2)uo«ik — +  2uoctW0(n-k)^ 

Jferz-  oo 


Using  an  identical  method  to  determine  a  solvability  condition  for  the  quantities 
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u ien  and  vic„,  we  also  find  that 


d°"V^  |  c 
Fif  2  +  ? 


1 _ X  (Q  +  2l/2)r  6 

3(o  +  2y2)  9xn/G^  3.T0en 


2  in  K\\/a  +  2y2 

3  •  y/Gox  a  +  2i/2 


^Ocn 


^Ocn^OcO 


+  (Nen!2  +  N,n|2)  ~  G0X(^  +  2?/2>0cn  “  ^yV  G()X(a  +  22/2) 

n=l 


-r  £  ( yj G0x(a  +  2y2)v0ck  —  °^~  —  ~  2i>o«t  W'oCn-t) 


(3.96) 

Finally,  to  close  the  system,  we  need  to  find  the  equations  which  determine  the 
harmonic  terms  U o„.  Considering  the  coefficients  of  En  in  (2.5a,  6),  we  find  that 


•  v  fr  .  dV°n  o 
tnK0Uon  +  — —  =  0, 


(3.10a) 


d2U0n 


.  (  ftSlN  2  dUly  1 

-in  [KlUO  +  -=■— J  Gon  -  ~zrv0n  7 ===== 


V  Gox  («  +  2?/2) 


(^Ocn^OcO) 


i  c?  r  °° 

-—=========—  53  (^0ckV0e(n-t)  +  V0,it>0»(n-*))  • 

2VGox(a  +  2y2)^  L  v 

(3.106) 

We  now  have  the  complete  specification  of  our  problem  to  determine  the  nonlinear 

A 

vortices.  The  steady  vortex  term  FoeO  is  given  by  (3.6),  the  mean  flow  Won  by  (3.7),  the 
wavy  mode  quantities  uoen>  vo,n  by  (3.9)  and  the  harmonic  terms  17 On  by  (3.10).  This 
system  comprises  an  eigenproblem  which  may  be  simplified  by  making  the  following 
transformations.  We  define  the  functions 


0l(®)  = 


6  x  {o.  +  2y2V 


ULV'  3(a  +  2j/2)  3  9xv/^X  ’ 

transform  the  independent  co-ordinate  £  to  £i  by  defining 

,  ,,92 f 
£l  =  €  +  — , 


«('>  =  v/C0X(3g  +  2”).  (3.11a, 6) 


(3.11c) 


write 


K{  =  Ki 


/«72^\/Gox 

<71  (a  +  2t/2)  2 


(3. lid) 
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and  finally  define  nondimensional  quantities  according  to 


v  V0e0  =  VG(-gi)'3A,  vocn  =  (-gi)' B0cn,  vo.n  =  {-gi)*Bo.n, 

WO*  (-<7l)’„  A  ^ 

Won  =  {-91) 3  Con,  "tt-  =  — : — Don,  n  = 


and 


Ki  = 


dr, 
Va  +  2  y2 

VGoX 


92 

(-9l)~  •  K\. 


(~9i)(a  +  2y2)  ’ 


(3.12) 


Inspection  of  the  resulting  equations  reveals  that  necessarily  we  must  have 


Bqc  (2j  - 1)  =  Bq,  (2j )  =  Cq(2j)  =  A)(2;-l)  =  0, 


for  all  integers  j,  and  then  if  we  define 


Boe(2j)  =  B2j,  -Go, (2,-1)  =  -82;- 1,  C'o(2,_i )  =  Cj,  D0(2j)  =  Dj,  (3.13 a) 

note  that  by  the  conjugacy  relationships  implied  by  (3.2)  we  have 


(3.13b) 


and  define  the  quantities  Bq  —  Dq  =  0  for  completeness,  we  obtain  the  system 
d2A 


dr,2 


—  tjA  —A3  + 


IB.P) 

r=l 


1  °° 

“h  53  [2-02,'-  l^l-i  ~  3B2jD„j  +  c.c.] , 


(3.14a) 


d2B2j-i 


2  i 


drj2 - VD2j-l  ~  j(2 j  -  l)(n?7  +  Ki)B2i-i  =  B2j-iA2 

+  ^Y1(E  |B.|2)  ~^§-C!A-\  E  (U32T.1Di-,  +  2B2,Cl.,), 


d2B2j  _  2i 


(3.14b) 


-  r,B2j  -  —(2j)(£lr)  +  Ki)B2j  =  52i  A* 


dr,2  '~"1  3 


+  ^  (E  \Br?j  -y/e DjA-^  E  (3B2rDj-r  —  252r_iC#_r-j:i) , 

r=l  r=-  00 

(3.14c) 
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(3.14  d) 


fin 

-i(2j-  l)(ClV  +  Kl)Cj 


\/6  (  Ad2B2j-i  „  d2A\ 
2  T  dr)  2  B2*-'dr?) 


•is(? 


and 


d2D< 


dr/2 
1  d2 


r  -  2ij(Slr)  +  K\)Dj  = 


3  dr)2  L2  _ 


1  “ 

r  (#2r#2;-2r  +  B2r_il?2i-2r+l)  +  \Z6B2jA 


(3.14c) 


As  yet  we  have  not  addressed  the  problem  of  specifying  the  boundary  conditions 
pertinent  to  the  system  (3.14).  Firstly,  we  anticipate  that  the  wavy  vortex  be  con¬ 
strained  to  lie  completely  within  the  shear  layer  Ila  and  this  requires  that  Bj,  Cj  and 
Dj  decay  exponentially  as  |?7|  — ►  oo.  Secondly,  the  steady  vortex  term  A  must  match 
with  the  core  region  solution  for  the  strongly  nonlinear  vortex  within  I.  This,  as  in 
HL,  implies  that  A  ~  (— rj)^  as  t)  — >  — oo.  Finally,  all  the  vortex  activity  is  confined 
within  the  core  J  and  the  shear  layer  Ila  so  that  in  order  that  in  Ilia  the  mean  flow 
is  governed  by  the  usual  boundary  layer  equations  we  impose  the  requirement  A  — *  0 
as  r)  — >  oo;  indeed  by  (3.14a)  we  have  A  oc  Ai(r))  in  this  limit  where  Ai  is  the  Airy 
function,  see  Abramowitz  &  Stegun  (1964). 

To  conclude  this  section,  we  note  that  the  nonlinear  wavy  modes  within  the  shear 
layer  Ila  are  determined  by  the  functions  Bj,  Cj,  Dj  which  are  defined  in  (3.12)  and 

(3.13) .  These  terms,  together  with  the  steady  vortex  quantity  A,  satisfy  equations 

(3.14)  with  associated  boundary  conditions 


and 


A  oc  Ai(r))  as  r)  — >  oo, 

(3.15a) 

A  ~  (~v)*  as  r) — >—oo, 

(3.156) 

Bj,Cj,Dj — >0  as  |?7| — >  oo. 

(3.15c) 

In  general,  the  solution  of  this  system  is  highly  complicated  numerical  task  and 
in  the  following  section  we  address  this  problem  in  a  weakly  nonlinear  limit  which 
permits  additional  analytical  progress. 
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§4.  Weakly  nonlinear  theory. 


We  recall  that  in  HS  the  infinitesimally  small  wavy  mode  problem  was  governed 
by  the  system  of  equations  (3.14)  for  the  functions  A,  Bi  and  C\  when  all  the  nonlin¬ 
ear  terms  (except  those  involving  the  function  A)  are  omitted.  HS  used  a  numerical 
scheme  to  search  for  eigenvalues  of  this  linearised  problem  but  located  only  one  such 
pair.  However,  HS  suggested  that  further  eigenvalue  pairs  with  K\  and  Cl  real  (and  so 
corresponding  to  neutrally  stable  wavy  modes)  may  be  possible  and  Bassom  &  Sed- 
dougui  (1989)  found  eight  such  pairs  in  all.  Further,  these  latter  authors  also  demon¬ 
strated  that,  plausibly,  infinitely  many  pairs  exist  and  they  derived  the  asymptotic 
structure  of  the  neutrally  stable  modes  in  a  high-wavenumber,  high-frequency  limit. 
The  eight  eigenpairs  identified  by  Bassom  &  Seddougui  (hereafter  BS)  are  summarised 
in  Table  1. 

Our  analysis  of  the  full  equations  (3.14)  is  concerned  with  the  problem  in  which 
the  wavy  mode  is  of  amplitude  0(h),  h  - Cl,  relative  to  the  scalings  implied  in  (3.2). 
In  this  case,  in  keeping  with  the  usual  weakly  nonlinear  approach,  we  anticipate  from 
(3.14)  that  the  perturbation  quantities  are  changed  by  0(h2)  from  their  linear  values, 
and  we  write 

A  =Ao  +  h2Ai  +  ...,  B\  =  hB\Q  +  h3Bu  +  . . . ,  i?2  =  h2B  20  +  . . . , 

C\  -hC\o  +  h3Cn  4-  . . . ,  C2  =  h2C20  +  •  •  • ,  D\  =  h2 D 10  4-  . . . ,  (4.1) 

Cl  =  (Cl)L  +  h2Clp  +  . . . ,  K\  =  (K\)l  +  h2Kp . 

The  remaining  terms  in  (3.14)  are  sufficiently  small  so  as  to  become  negligible  in  the 
ensuing  analysis. 

Substituting  (4.1)  in  (3.14)  and  (3.15)  gives,  at  leading  orders  in  h,  the  linear 
eigenvalue  problem  described  by  HS.  In  this  paper,  together  with  BS,  figures  of  the 
linear  eigenfunctions  B 10  and  C10  appropriate  to  some  of  the  eigenpairs  in  Table  1  are 
shown.  We  notice  from  (3.14a),  (3.15a,  6)  that  the  steady  vortex  profile  Ao  satisfies 

<^Aq  .  .3 

~drj^  ~  VA0  ~  j4°’ 

j4o — >A Ai(ii)  as  v — >00,  Ao  ~  (— 77)’  as  rj — >  —00. 

Further,  the  results  of  Hastings  &  McLeod  (1980)  demonstrate  that  A  =  \/2  and  then 
the  function  Ao  may  be  evaluated  by  using  the  asymptotic  condition  as  77  — >  00  and 
integrating  in  the  direction  of  decreasing  77  using  a  Runge-Kutta  method,  see  HS. 
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At  next  order  in  h  when  (4.1)  is  inserted  in  (3.14),  we  obtain 

/  ^2  \  12 
V  dn 2  —  *0  Al~  3A0A1  =  3  ^-ol-Slol2  +  {B\qC{q  +  c.c .)  , 


drj t 

+  ai{'n)Bn  +  a2(v)Cn  =  AO?)  +  KP£2(v)  +  £3(7?), 

d2Cn  ,  ,  ^  ,  ,  sr,  r>  f-  ,  \/6Ao£i(?7) 

-j~2 — ^03(77)^11  +  a4(7?)Cn  =  Sip  |^7?Cio  4 - - - 


+  I(V 


iClO  + 


VQAoC2(v) 


+  ^-CsW  +  AM, 


(4.2  a) 
(4.26) 

(4.2c) 


where 


ai  =  —  V  1  + 


2  in  1 


2i(^l)i  ,2  2\/6 

- - - A0,  a2  =  — A0, 


(4.3a,  6) 


a3 


—  —  ~ r -  +  {Ki)l')  Ao,  04  =  — i  (0^7?  +  (I^i)l)  +  2Aq.  (4.3c, d) 


Here,  we  find  that  C\, . . . ,  £4  are  defined  by 

2i  2i 

£l  =  —  vBio,  £2  =  y-BlO, 


(4.4a, 6) 


£3  =  2AoAii?io  + 


-Biol^iol2  2\/6 


3  C10A1  -  B*10Dio  -  -J?2oCio> 


(4.4c) 
■HloAol^iol2 


£4  =  -  V^AiAq^io  -  2A0A1C10  +  -y-  [yi'L’n  +  (■K’l)  ^ )  A1B10 - —  ~^= 

-  ^  (5ioCio  +  2£i0C10)  -  i  {nLv  +  (k 1)  )  5i0520 

O  Xr 

,  DioAoBloVe 
-  —AoB2oC10  +  - - - , 

(4.4d) 

where,  from  HS,  we  know  that  the  functions  I?io  and  C3.0  satisfy  the  homogeneous 
forms  of  equations  (4.26,  c).  Additionally,  from  (3.14)  we  obtain 

d  fo°-  ~  VB20  ~  %  (&lV  +  (ki)  )  B20  -  AqB20  =  I^iqCio  ~  V6D10A0,  (4.5 a) 
a?7z  3  '  N  l  3 


and 


-  2 i  (6*17  +  (kx)  J  D10  =  +  V6B20A0]. .  (4.0 b) 
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To  ensure  the  disturbance  is  confined  within  the  shear  layer  we  demand  A\ ,  Du, 
Cn,  B20  and  Dio  — >  0  as  |?7|  — >  00. 

The  homogeneous  forms  of  (4.26,  c)  are  solved  in  HS.  Then,  as  is  the  normal 
procedure  within  a  weakly  nonlinear  analysis,  (4.26,  c)  only  have  a  suitable  solution 
if  a  certain  compatibility  condition  is  satisfied.  This  condition  leads  to  the  speci¬ 
fication  of  the  correction  terms  Kp  and  within  the  wavenumber  and  frequency 
expansions  detailed  in  (4.1).  To  derive  our  compatibility  requirement  we  first  need 
to  consider  the  system  adjoint  to  the  homogeneous  forms  of  (4.26,  c).  The  adjoint 
functions  (G  1(7]),  G 2(1]))  satisfy  the  coupled  equations 

d?Gi 

+  ai(r?)Gi  +  a2(v)G2  =  0,  (4.6a) 

3T  +  a2{v)Gi  +  04(77)^2  =  0,  (4.66) 

with  associated  boundary  conditions  of  G 1,  G2  — >  0  as  I77I  — >  00.  Multiplying  (4.26) 
by  Gi(tj),  (4.2 c)  by  (*2(77),  adding  the  resulting  equations  and  integrating  by  parts, 
reveals  that  for  a  satisfactory  solution  of  (4.26,  c)  to  exist,  we  require 

/°°  {G1N1  +  G2N2)dr1  =  0,  (4.7) 

-  00 

where  N\  and  N2  are  the  right  hand  sides  of  (4.26,  c)  respectively. 

On  substituting  (4.4)  into  (4.7)  we  eventually  obtain  a  relationship  between  the 
wavenumber  correction  Kp  and  the  frequency  correction  flp  of  the  type 

1  +  ziflr  +  z2Kp  =  0,  (4.8) 

where  z\,  z2  are  complex  numbers  expressible  as  quotients  of  integrals.  In  particular, 
z\  =  /1//3  and  z2  =  I2//3  where 

h  =  /”  [,CiGl+  (»)ClO  +  G2  in,  (4.9a) 

h  =  j'°  [.C2G1  +  (iCio  +  V^°£2)  <?2  in,  (4.96) 

and 

h  =  I  ”  [ -CsGi  +  (  (-  c, |  G2  in-  («°) 
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Before  presenting  results  of  the  computations,  we  observe  that  we  may  easily 
generalise  our  analysis  to  allow  the  wavy  perturbation  to  evolve  on  a  slow  lengthscale 
so  that  the  growth  or  decay  of  a  non-neutral  disturbance  may  be  monitored.  Formally, 
if  we  consider  the  neighbourhood  of  a  point  x  =  xn  (at  which  an  infinitesimal  wavy 
vortex  is  neutrally  stable)  and  allow  a  constant  frequency  disturbance  to  develop  on 
the  lengthscale  X,  where  X  =  e~*h~2(x  —  xn ),  h  <C  1,  then  the  result  of  repeating 
the  previous  analysis  for  a  disturbance  of  amplitude  hA(X)  relative  to  the  scalings  in 
(3.2)  is  the  evolution  equation 


(1A 

dX 


J2,a  A- -\A? A. 

Z2  22 


(4.10) 


We  may  then  easily  obtain  the  equation  for  the  development  of  the  wavy  disturbance 
amplitude 

(|i|2)  =  2Re  (^ji)  n„|i|2  +  2Re  (=^)  |i|4.  (4.11) 

To  effect  the  necessary  calculations  we  used  the  method  described  by  HS  to  ob¬ 
tain  the  linear  eigenfunctions  B 10  and  C\q.  To  fix  a  normalisation  for  the  problem,  we 
elected  to  scale  these  eigenfunctions  such  that  max{|i?io|  :  |r?|  <  00}  =  1.  The  tech¬ 
nique  of  HS  is  easily  adapted  to  solve  the  adjoint  system  (4.6)  and  essentially  consists 
of  the  following.  Equations  (4.6)  were  written  as  a  set  of  four  first  order  differential 
equations,  which  was  solved  by  using  a  fourth  order  Runge-Kutta  scheme.  Since  (4.6) 
is  the  system  adjoint  to  the  problem  specified  in  HS,  the  two  problems  have  the  same 
e:genvalucs  and  the  integration  procedures  for  (4.6)  were  started  at  77  =  —00  and  at 
77  —  00  and  both  continued  to  77  =  0.  The  detailed  asymptotic  conditions  on  G\  and 
G 2  for  large  |t7|  follow  in  a  similar  manner  to  those  given  in  HS  and  so  yield  the  ini¬ 
tial  conditions  for  the  integrations.  The  two  integration  procedures  from  ±00  were 
matched  across  77  =  0,  and  hence  the  adjoint  functions  G\,  G2  determined. 

In  order  to  obtain  the  integrands  in  (4.9)  we  needed  to  find  the  functions  Ai,  i?20 
and  Z>io  in  addition  to  those  already  computed.  The  steady  vortex  quantity  Ai,  which 
satisfies  |Ai|  — >  0  as  [77!  — >  00  was  evaluated  using  a  finite  differencing  procedure  on 
(4.2a).  We  observe  from  (4.5)  that  the  remaining  wavy  vortex  terms  D 10  and  B 20  are 
governed  by  coupled  equations  and  these  functions  were  found  also  by  employing  finite 
difference  methods.  Once  these  flow  quantities  were  all  determined,  the  necessary 
integrations  to  evaluate  (4.9)  were  effected  by  use  of  trapezoidal  integrations  and 
the  results  so  obtained  were  checked  to  ensure  that  they  were  independent  of  tin. 
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step  lengths  used  within  the  Runge-Kutta  and  finite  differencing  procedures  and  also 
independent  of  the  (large)  values  of  |r?j  at  which  the  various  asymptotic  behaviours  as 
| Ty |  — >  oo  were  assumed  to  be  valid. 

The  results  of  these  computations  and  their  implications  are  considered  in  the 
following  section. 


§5.  Results  and  discussion. 


The  numerical  scheme  outlined  above  was  used  to  determine  the  coefficients  in 
the  amplitude  equation  (4.11).  This  equation  demonstrates  that  the  weak  nonlinearity 
of  the  problem  allows  the  existence  of  a  threshold  equilibrium  amplitude  At  given  by 


Ret-.WA  „ 

Ref-izJ1)  J  r 


(5.1) 


The  relevant  values  of  the  quantities  contained  within  (5.1)  are  summarised  in  Table 
2,  where  we  list  those  corresponding  to  the  first  five  linearly  neutrally  stable  modes 
given  in  Table  1.  Eigenfunctions  for  the  problem  (4.2)-(4.5)  in  the  cases  ClL  =  0.372 
and  Cl i  =  0.659  are  shown  in  Figures  2  and  3  respectively. 

Inspection  of  the  complex-valued  eigenfunctions  reveals  that  the  imaginary  parts 
are  pictorially  similar  to  the  real  parts  and  so  we  elect  to  concentrate  upon  the  latter 
in  presenting  the  eigenfunctions  of  the  problem.  Observation  of  the  eigensolutions  for 
the  first  five  linearly  neutrally  stable  modes  shows  that  the  eigenfunctions  G i,  C?2,  B20, 
£>10  have  the  same  trends  as  the  linear  eigensolutions  Bio,  Cio  described  by  Bassom 
&  Seddougui  (1989).  Thus,  the  eigenfunctions  are  oscillatory  in  nature  and  centred 
about  a  negative  r)  position.  Moreover,  this  position  is  close  to  rj  =  —(K\)LlClL.  We 
also  note  that  the  number  of  oscillations  of  the  complex  eigenfunctions  increases  for 
each  eigenmode  with  the  eigenfunctions  corresponding  to  the  choice  = 

(0.690,0.372)  having  the  least  number  of  oscillations. 

Since  the  coefficient  of  the  nonlinear  term  in  (4.11)  is  negative  for  each  of  the 
cases  considered  in  Table  2,  interest  is  concentrated  on  supercritical  disturbances  when 
np  <  0.  Then,  on  the  basis  of  linear  theory,  the  wavy  modes  of  Table  2  are  unstable, 
but  the  nonlinearity  in  (4.11)  ensures  that  the  equilibrium  amplitude  At  is  stable  in 

A  A  A  A 

the  sense  that  for  A  <  Ae  the  disturbance  will  tend  to  grow  whereas  for  A  >  At  it 
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will  decay  until  A  =  Ac.  Of  course,  liere  we  have  only  discussed  the  results  for  the 
five  modes  of  lowest  frequencies  and  from  the  work  of  Bassom  &  Seddougui  (1989)  we 
know  that  it  is  probable  that  there  are  infinitely  many  more  modes.  However,  it  is 
the  case  that  in  the  development  of  the  flow  over  a  curved  wall  the  lower  frequency 
vortices  are  likely  to  be  the  more  important:  an  aspect  which  we  consider  in  more 
detail  later. 

The  procedure  described  above  sought  a  description  of  the  evolution  of  a  wavy 
mode  in  terms  of  the  perturbation  ftp  of  the  non-dimensional  frequency  away  from  the 
linear  neutral  values  (ft)  .  We  remark  that  in  order  to  translate  our  findings  for  the 

h 

case  of  any  particular  curved  wall  of  prescribed  curvature  x(a:)>  we  may  only  determine 
the  dimensional  wavenumber  I(\  and  frequency  ft  once  the  scaling  quantities  in  (3.11) 
and  (3.12)  have  been  calculated.  To  do  this  requires  an  application  of  the  computation 
described  by  HL  for  the  particular  x(x )•  In  addition  to  their  account  of  this  numerical 
calculation,  HL  investigated  two  asymptotic  limits.  Firstly,  solutions  were  determined 
for  their  free-boundary  problem  for  x  close  to  the  linear  neutral  position  x  =  x.y  and 
secondly,  for  x  far  downstream  of  this  point.  In  the  latter  case  it  was  assumed  that 
the  curvature  x  increases  at  least  as  quickly  as  x »  for  x  »  1. 

Rather  than  attempt  to  pursue  the  HL  calculation  for  any  specific  boundary  layer 
flows  here,  it  was  felt  that  it  would  be  more  instructive  to  draw  general  conclusions 
concerning  nonlinear  wavy  vortices  which  would  hold  for  large  classes  of  flows.  This  is 
largely  motivated  by  the  fact  that  all  the  experimental  papers  relating  to  wavy  vortices 
within  boundary  layer  flows  over  curved  walls  contain  insufficient  detail  to  enable  the 

A  A 

dimensionless  values  of  K\  and  ft  to  be  computed.  Hence  any  results  we  may  draw 
concerning  specific  flows  could  not  be  easily  related  to  practical  observations. 

We  now  turn  to  consider  the  behaviour  of  the  frequency  of  the  neutrally  stable 
nonlinear  modes  in  each  of  the  asymptotic  limits  considered  by  HL.  As  x  — *  xt, 
where  we  recall  the  definition  that  x. -is  the  location  where  the  steady  Gortler  vortex 
is  linearly  neutrally  stable,  we  know  that  the  two  shear  layers  merge.  Then,  following 
the  notation  of  HS,  if  we  denote  by  ftT  and  f lB  the  (dimensional)  frequencies  of  the 
neutral  modes  in  the  upper  and  lower  shear  layers  respectively,  we  have  from  (3.12) 
that  as  x  — >  x , ,  ftT  — >  fts .  We  now  suppose  that  for  large  xf  /v  xM ,  M:  >  j.  This 
choice  of  curvature  is  again  chosen  to  illustrate  the  type  of  beliayiour  we  may  observe 
in  a  practical  situation  rather  than  as  representing  any  one  particular  boundary  layer 
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flow.  In  this  case,  the  behaviour  of  the  functions  a,  yi,  y2,  b  are  as  given  by  HL  and 
by  routine  application  of  these  forms,  we  may  easily  show,  as  in  HS,  that  for  x  »  1, 

nT  ~  ^-mgI x2*1-1,  ~  in.  (5.2) 

y  o 

Consequently,  as  x  increases,  the  frequency  of  the  neutral  wavy  mode  confined  to  the 
upper  shear  layer  increases  whereas  that  of  the  mode  in  the  lower  layer  tends  to  a 
constant  value.  Thus  the  modes  can  be  categorised  as  being  of  high  and  low  frequency 
respectively. 

We  remark  that  the  instability  process  detailed  in  section  3  for  the  upper  layer 
is  also  perfectly  possible  within  the  lower  layer  lib,  and  if  exponentially  small  terms 
may  be  neglected,  the  two  shear  layers  are  completely  independent  for  the  unstable 
travelling  waves  decay  exponentially  away  from  the  centres  of  Ila ,  b.  Theoretically  at 
least,  our  analysis  could  be  adapted  to  prove  that  an  0(1)  disturbance  in  one  shear 
layer  provokes  an  exponentially  small  response  in  the  other. 

In  a  practical  situation,  it  is  often  found  that  a  wavy  mode  of  a  predetermined 
frequency  is  imposed  upon  the  flow  and  the  development  of  the  disturbance  monitored 
as  it  moves  downstream.  Then  we  anticipate  that  the  two  shear  layers  will  breakdown 
at  different  ^-locations.  As  discussed  by  HS,  since  the  downstream  velocity  component 
of  the  the  basic  state  is  larger  within  the  upper  layer  Ila ,  it  is  expected  that  this  layer 
is  the  first  to  become  unstable.  Further,  as  the  wavenumber  K(x )  relating  to  a  fixed 
frequency  perturbation  will  be  different  in  the  two  shear  layers,  the  modes  in  the  layers 
move  with  different  wavespeeds.  Close  to  the  linear  neutral  position,  the  shear  layers 
meet  and  so  if  breakdown  occurs  close  to  this  point  then  the  structures  within  the  two 
layers  will  be  very  similar.  Far  downstream,  the  layer  Ila  moves  into  the  freestream 
and  the  downstream  velocity  component  tends  to  the  freestream  speed.  The  lower  layer 
sinks  towards  the  bounding  wall  so  the  fluid  velocity  there  approaches  zero.  Thence 
if  the  stationary  vortices  develop  over  a  sufficiently  lengthy  interval  before  breakdown 
occurs,  the  mode  in  the  upper  layer  has  speed  close  to  that  of  the  freestream,  whilst 
the  one  in  lib  has  a  very  small  speed. 

We  now  consider  the  implications  of  our  analysis  for  an  experimental  setting. 
We  note  from  (5.2)  that  for  the  case  of  a  mode  of  prescribed  dimensional  frequency  0 
imposed  upon  the  flow  upstream  of  the  linear  neutral  point,  as  we  move  far  downstream 
(for  the  case  of  x  ~  xM)  the  non-dimensional  frequency  Cl  for  vortices  confined  to 
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the  upper  layer  Ila  tends  towards  zero.  Our  results  obtained  here  demonstrate,  as 
previously  mentioned,  the  existence  of  a  stable,  threshold  equilibrium  amplitude,  at 
least  for  the  weakly  nonlinear  modes  close  (in  parameter  space)  to  one  of  the  first 
five  linearly  neutral  disturbances.  As  the  perturbation  propagates  downstream,  the 
non-dimensional  frequency  Cl  of  modes  trapped  within  the  upper  shear  layer  decreases 
and  our  analysis  indicates  that  there  are  no  unstable  non-zero  equilibrium  amplitudes 
within  this  frequency  regime.  Hence,  we  would  anticipate  that  the  low-frequency  (low 
Cl)  modes  investigated  in  this  paper  are  likely  to  be  the  most  important  in  a  practical 
experiment  and  that  finite  amplitude  wavy  vortices  are  likely  to  be  observable  in 
practice. 

In  this  paper  we  have  chosen  not  to  pursue  an  extensive  solution  of  (3.14)  away 
from  the  domain  where  the  weakly  nonlinear  theory  is  applicable.  In  such  a  calcula¬ 
tion,  all  the  unknowns  in  (3.14)  are  potentially  important  and  would  necessarily  have 
to  be  accounted  for  in  a  fully  rational  way.  Preliminary  studies  of  the  computation  of 
the  full  system  (3.14)  and  (3.15)  suggest  that  this  task  is  far  from  easy.  It  would  be 
interesting  to  make  a  careful  examination  of  these  equations  to  determine  the  char¬ 
acteristics  of  strongly  nonlinear  wavy  modes.  However,  here  we  have  elected  to  deal 
only  with  the  weakly  nonlinear  approach.  The  results  of  this  analysis  suggest  that 
for  numerical  computations  based  on  the  full  version  of  (3.14)  it  can  be  speculated 
that  for  each  non-dimensional  low  frequency  Cl  there  is  a  multitude  of  possible  vortex 
states,  each  stable  and  with  a  different  wavelength.  Equivalently,  for  a  wavy  vortex 
of  prescribed  streamwise  wavenumber,  there  could  conceivably  be  a  family  of  possi¬ 
ble  non-dimensional  frequencies  for  this  disturbance.  A  thorough  description  of  the 
possible  solutions  of  (3.14)  is  eagerly  awaited. 

Papers  by,  among  others,  Bippes  &  Gortler  (1972),  Wortmann  (1969)  and  Bippes 
(1978)  describe  some  experimental  investigations  of  the  secondary  instability  of  Gortler 
vortices.  Recently,  work  by  Kohama  (1987)  and  Peerhossaini  &  Wesfreid  (1988)  have 
given  more  details  concerning  this  secondary  breakdown  than  did  the  aforementioned 
authors.  Kohama  (1987)  performed  experiments  on  a  NASA  laminar  flow  wing.  He 
considered  the  breakdown  only  in  the  upper  parts  of  the  vortices  and  gave  measure¬ 
ments  of  the  speeds  of  the  wavy  modes  at  various  downstream  locations.  Peerhossaini 
&  Wesfreid  (1988)  examined  a  boundary  layer  over  a  concave  section  of  a  curved  chan¬ 
nel.  They  reported  that  when  the  secondary  instability  first  manifested  itself,  it  did 
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so  at  the  top  of  the  steady  Gortler  vortices  and  that  further  downstream  another  sec¬ 
ondary  mode  was  observed  near  the  wall.  Both  instabilities  had  wavy  boundaries  and 
so,  like  HS,  we  conclude  that  these  modes  are  formed  as  a  result  of  the  mechanism 
considered  in  this  paper.  This  experimental  work  of  Peerhossaini  &  Wesfreid  (1988) 
lends  credence  to  our  conclusion  that  according  to  weakly  nonlinear  theory  the  finite 
amplitude  wavy  vortices  are  stable  and  hence  expected  to  be  obtainable  in  practice. 
Interestingly,  Peerhossaini  &  Wesfreid  also  commented  that  as  the  disturbance  devel¬ 
oped  downstream  no  coherent  frequency  of  the  mode  was  apparent.  Possibly,  in  this 
case,  the  weakly  nonlinear  approach  is  no  longer  tenable  and  the  full  system  (3.14)  is 
appropriate;  we  then  have  the  possibility  of  the  existence  of  a  whole  family  of  modes 
with  different  frequencies  as  discussed  above.  This  may  explain  the  experimental  ob¬ 
servation  of  Peerhossaini  &  Wesfreid  (1988). 

To  conclude,  we  observe  that  we  have  shown  that  according  to  our  weakly  non¬ 
linear  theory,  for  non-dimensional  frequencies  close  to  the  linear  neutral  values,  finite 
amplitude,  stable  equilibrium  wavy  modes  exist.  Our  results  are  consistent  with  the 
available  experimental  observations;  however,  the  lack  of  detail  given  in  the  accounts  of 
this  practical  -work  means  that  precise  comparisons  between  those  studies  and  our  find¬ 
ings  is  rendered  near  impossible.  Of  interest  would  be  the  consideration  of  the  strongly 
nonlinear  system  (3.14)  although  we  ourselves  are  not  planning  to  pursue  this  aspect 
further.  However,  we  do  feel  that  this  investigation  merits  attention  at  some  stage. 
We  finally  note  that  other  instability  processes  could  give  rise  to  three-dimensionality 
and  time-dependence  in  the  latter  stages  of  vortex  development;  these  would  most 
likely  be  Rayleigh  instabilities  associated  with  the  spanwise  locally  inflexional  velocity 
profiles,  or  Tollmien-  Schlichting  waves,  both  of  which  need  to  be  thoroughly  studied 
to  enable  a  more  complete  picture  of  the  transition  process  to  be  constructed. 
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Table  1.  Values  of  K\  and  Cl  for  which  infinitesimally  small  wavy 
vortex  modes  are  neutrally  stable,  from  Bassom  &  Seddougui  (1989). 


mm 

Kfl 

-2izi/z2 

-2i/z2 

0.6895 

0.3715 

IS&fcS&ill 

(-0.187,0.345) 

3.89 

2.900 

0.6587 

(-1.019,  -6.089  xl0“2) 

2.033 

4.156 

0.7416 

(-1.044,-1.720) 

2.325 

5.435 

0.7950 

(-7.06,10.09) 

(-0.121,-3.997) 

7.652 

7.527 

0.9976 

(-4.57,11.86) 

(-1.335,-4.461) 

1.849 

Table  2.  Quantities  defined  in  (4.11)  and  (5.1)  appropriate 
to  the  weakly  nonlinear  analysis  of  section  4. 
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FIGURE  CAPTIONS 


Figure  1.  The  different  regions  beyond  the  downstream  position  of  neutral  stabil¬ 
ity. 


Figure  2.  Eigenfunctions  for  the  weakly  nonlinear  problem  of  §4  as  defined  by 
(4.1)-(4.4)  when  (ICi)t  =  0.690,  ClL  =  0.372.  (a)  Re{Bwy  ( b )  Re{C10)\  ( c )  Re(Gi); 
(d)  Re(G2y,  (e)  AV)  (/)  Re(B20 )  and  ( g )  Re(D10). 


Figure  3.  Eigenfunctions  for  the  weakly  nonlinear  problem  of  §4  as  defined  by 
(4.1)-(4.4)  when  (lti)L  =  2.900,  =  0.659.  (a)  Re{B10y  (6)  Re{C10)-,  (c)  Re{G{)\ 

(d)  Re(G2)-,  (e)  AV,  (/)  Re(B2 0)  and  ( g )  Re(Dio). 
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Figure  3f 
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